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^^ ■ We present a six-dimensional potential energy surface for the (H2)2 dinier based 

>. 

O ■ on coupled-cluster electronic structure calculations employing large atom-centered 

(N 



Gaussian basis sets and a small set of midbond functions at the dimer's center 
of mass. The surface is intended to describe accurately the bound and quasibound 



r~| . states of the dimers (H2)2, (02)2; and H2-D2 that correlate with H2 or D2 monomers 

Oh: 

in the rovibrational levels {v,j) = (0,0), (0, 2), (1, 0), and (1, 2). We employ a close- 

D ' coupled approach to compute the energies of these bound and quasibound dimer 

^ ' states using our potential energy surface, and compare the computed energies for 

. ^ , infrared and Raman transitions involving these states with experimentally measured 

^~>' transition energies. We use four of the experimentally measured dimer transition 

^ ; 

I *" <■ energies to make two empirical adjustments to the ab initio potential energy surface; 

the adjusted surface gives computed transition energies for 56 experimentally ob- 
CN served transitions that agree with experiment to within 0.036 cm . For 29 of the 56 

cn ■ transitions, the agreement between the computed and measured transition energies 

is within the quoted experimental uncertainty. Finally, we use our potential energy 

r^ , surface to predict the energies of another 34 not-yet-observed infrared and Raman 

O 

transitions for the three dimers. 
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I. INTRODUCTION 

The (H2)2 dimer has long been viewed as a prototypical bimolecular van der Waals dinier. 
Because the (H2)2 dinier is electronically simple, it has been the focus of a number of ab 



initio studies;-'2ii^i^i^i^>^'^ii2i>iiii2i>i^i^ii^i^>i^ however, because the H2-H2 van der Waals inter- 
action is quite weak,— ab initio calculations with accuracy much higher than the oft-quoted 
"chemical accuracy" of 1 kcal/mol must be employed to provide useful information about 
the H2-H2 potential energy surface. Two recent advances in ab initio methods have made 
it possible to compute the H2-H2 interaction with the required level of accuracy: (1) the 
development of hierarchical sequences of one-electron Gaussian basis sets for approximating 
molecular electronic wave functions,— sequences which systematically approach the complete 
one-electron basis set limit, and (2) the development of efficient methods for accounting for 
electron correlation effects in these wave functions by systematically approaching the many- 
electron basis set limit.— "^ 

Diep and Johnsonr^>^ took advantage of these two advances in ab initio methods to 
compute an accurate four- dimensional rigid rotor potential energy surface for the (H2)2 
dimer; the Diep- Johnson surface gives low temperature second virial coefficients and integral 
elastic scattering cross sections in reasonably good agreement with experiment. However, 
this potential energy surface does not depend explicitly on the covalent bond lengths of 
the two H2 monomers, and thus is only able to describe the interaction between two H2 
molecules in their v = vibrational ground states. 

More recently, Boothroyd et al}^ have compiled a large database of energies for the H4 
system, based largely on multireference configuration interaction ab initio calculations, and 
have fit a global six-dimensional H4 potential energy surface to these energies. However, 
in this database, (H2)2 dimer configurations representative of the van der Waals well are 
assigned energies that come not from ab initio calculations, but rather from a empirically 
modified rigid rotor potential energy surface. Recent theoretical studies of low-energy in- 
elastic H2-H2 collisions that use this potential energy sniia.c^'^^^^ yield computed energy 
transfer rate coefficients in rather poor agreement with experiment. 

The (H2)2 dimer has also been the focus of several experimental investigations, beginning 
with the pioneering work of Watanabe and Welsh-- that demonstrated the dimer's existence 
through observation of its infrared (IR) absorption spectrum in the H2 f = 1 <— vibra- 
tional fundamental band. Later experimental studiea^^i^l recorded at high resolution the IR 
absorption spectra of the (H2)2 dimer (and several of its isotopomers) in the f = 1 ^ fun- 
damental band and f = 2 ^ first overtone band of the corresponding monomers. The high 
resolution IR absorption spectra of (H2)2 in the H2 fundamental and overtone regions, and 



the analogous isotopomer spectra, provide information about the vibrational dependence of 
the H2-H2 interaction. Complementary studies^^ of the far-IR absorption spectrum of the 
dimer provide information about the anisotropy of the potential energy surface in the region 
of the van der Waals well. 

Recently, the Raman spectrum of the (H2)2 dimer in the H2 fundamental region has also 
been observed. ^^ This spectrum provides information about the vibrational dependence of 
the H2-H2 interaction that is complementary to that provided by the high-resolution IR 
studies. Specifically, the vibrationally excited state of (H2)2 that is probed by the IR studies 
is one in which the vibrational excitation is delocalized across the two H2 monomers in 
an antisymmetric fashion, while in the Raman studies, the excited (112)2 state is one in 
which the vibrational excitation is delocalized symmetrically across the two monomers. A 
comparison of the IR and Raman spectra thus provides insight into the coupling between 
the two H2 vibrational modes in the (112)2 complex and into the dependence of the H2-H2 
potential energy surface on the two monomers' bond lengths. 

Equipped with this new information, we attempt here the construction of a six- 
dimensional II2-H2 potential energy surface that accurately describes the dimer's van der 
Waals well. We begin by computing ab initio H2-H2 interaction energies that are nearly 
converged with respect to both the one-electron and many-electron basis sets, and then 
construct a smooth potential energy surface from these computed interaction energies. We 
then make two small empirical adjustments to the surface; these adjustments soften slightly 
the surface's short-range repulsive wall, and increase slightly the strength of the surface's 
anisotropic term that couples the rotational degrees of freedom of the two monomers. The 
empirically adjusted surface gives IR and Raman transition energies for the {para-B.2)2, 
{ortho-D2)2i and para-Y{2-ortho-D2 dimers in good agreement with available experimental 
data.^'^s^ 

II. AB INITIO COMPUTATIONS 

A. Functional form of the H2 H2 interaction 

We consider a space-fixed coordinate system (x, y, z) in which one H2 molecule (denoted 
molecule 1) has its center of mass at the origin and the other H2 molecule (denoted molecule 



2) has its center of mass on the positive z axis. The orientation of molecule i is specified by 
its spherical polar and azimuthal angles (6*4, 0j). We let R represent the distance between the 
molecules' centers of mass, and let r, represent the bond length of molecule i. The H2-H2 
potential energy surface can then be expanded in terms of coupled spherical harmonics:— 

y=Y. ^h,hAR, n, r2)Gi,,i,,L{9,, 02, 0) (1) 

h,l2,L 

where = 02 — 0i, the summation indices /i, I2, and L are non-negative integers that must 
satisfy 

h + h + L = even integer, (2) 

and the homonuclear symmetry of the two H2 monomers dictates that h and I2 are also 
both even. The angular functions Gi^^i^^l have the form 



Gi^,i^,L = y^ ^C{li,m,l2,-m]L,0)Yi^^rn{0i,(pi)yh'm{0i,(p2) (3) 
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where C is a Clebsch-Gordan coefficient and Yi^m is a spherical harmonic normalized so that 
^,m(0,0) = (5m,o\/(2/ + l)/47r. (We use the Condon-Shortley phase convention for Yi^m-) 
The appearance of the Clebsch-Gordan coefficient C in Eq. ([3]) means that h, I2, and L 
must satisfy the angular momentum triangle rule. 

The functions Gi^^i^^l constitute a complete, orthogonal basis set for functions of the three 
angular coordinates {61,62, (p). For fixed R, ri, and r2, the coefficient Ai^^i^^i^R^r 1,1-2) can 
therefore be computed as 

K,i,,l{R. ri, r2) = ^^ j j Gi,,i,,l{6,, 62, ^)ViR, n, r2, 6,, 62, 0) d^i dS2 (4) 

where dSi = sin 6i d6i d0j. 

Earlier studies of the four- dimensional rigid-rotor H2-H2 potential energy surface^'^^'^^ 
show that the surface is dominated by four terms, with (/i, I2, L) = (0, 0, 0), (0, 2, 2), (2, 0, 2), 
and (2,2,4). In this work, we use numerical quadrature to compute the right-hand side of 
Eq. (jlj) for these four (/1,/2,-L) triples. Specifically, at fixed values of R, ri, and r2, we use 

to evaluate the integrals 



the 18-point spherical quadrature rule numbered 25.4.64 in Ref. 131 
over both dSi and d5'2 in Eq. (jl]). This requires us to compute the H2-H2 interaction energy 
K(-R, ri,r2, 6*1, 6*2, 0), using ab initio quantum chemical methods that we describe in the 
next subsection, at 12 sets of angles {6i,62,4>) when ri = r2 and at 19 sets of angles when 



Ti ^ r^- Symmetry relationships allow the rest of the 18^ = 324 interaction energies at fixed 
(i?, ri, r'2) to be determined from these ab initio calculations. 

The accuracy of the Ai^^i^^i^ coefficients computed in this fashion is limited by the fact 
that the quadrature rule we use fails to reproduce the orthogonality conditions 



y y G'ii,/2,L(^i, ^2, 0)G',i,/^,L'(^i, ^2, 0) d5i d52 = bi.^vM'.^^-^^'^^^ 



(5) 



when l\ + /'i > 6 or /2 + ^2 — 6- This means that the value of Ao,o,o obtained via quadrature 
also includes some contamination from v46,o,6 and v4o,6,6 (if these coefficients are nonzero in 
the ab initio potential energy surface), while ^2^2,4 is contaminated by (among other terms) 
^2,4,6 and A4 2,65 which describe the long-range electrostatic quadrupole-hexadecapole (QH) 
interaction between the two H2 molecules. 

To assess the magnitude of these erroneous contributions to the four Ai-^j^j coefficients of 
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to calculate 



interest, we used the more accurate 24-point spherical quadrature rule of Ref. 
the coefficients at (i?, ri,r2) = (4.5 oq, 1.4 oq, 1.7 oq), a repulsive (132)2 configuration where 
we expect the angular anisotropy of the potential energy surface to be relatively high, and 
where this contamination should thus be relatively severe. Table [T] compares the coefficients 
obtained using the two quadrature rules (based on ab initio interaction energies computed 
using the protocol outlined in Sec. Ill Bj) : the errors introduced at this {R, ri, r2) configuration 
by using 18-point quadrature appear to be quite small for the four terms that we include 
in our final potential energy surface. This table also gives the values for two additional 
coefficients in the coupled spherical harmonic expansion, ^2.2,0 and ^2.2,2? at this (1^2)2 
configuration, and shows that they are one to two orders of magnitude smaller than any of 
the four terms we retain in Eq. ([T]). This is in accord with previous studies^'^'^ of the 
four-dimensional rigid-rotor H2-H2 potential energy surface. 

B. CCSD(T) ab initio calculations 

We use Gaussian 03^ to compute the H2-H2 interaction energy, employing a coupled- 
cluster— >^ treatment of electron correlation that includes single and double excitations and 
a perturbative treatment of triple excitations,^'^ abbreviated CCSD(T). The CCSD(T) cal- 
culations are based on a restricted Hartree-Fock reference wave function; we have verified 
that such a reference does not exhibit a restricted — > unrestricted instability for the H2 bond 
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lengths considered here. We use the aug-cc-pVQZ basis set^i^ for the four hydrogen atoms, 
supplement this atom-centered basis set with a set of (3s3p2d) bond functions positioned at 
the dimer's center of mass, and employ the standard counterpoise correction.— The bond 



function exponents are taken from Ref. 
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We carry out these calculations at ri and r2 values of 1.1, 1.4, and 1.7 ag, and at 19 R 
values ranging from R = 4.25 Qq to 12.0 Oq, for a total of 1653 unique {R,ri,r2,0i,92,4>) 
(H2)2 configurations. We turn off automatic checking of the one-electron overlap matrix for 
near linear dependence and retain all 206 one-electron basis functions at every configuration; 
this eliminates possible discontinuities in the potential energy surface that could arise when 
some of these functions are dropped from the one-electron basis set. The Gaussian 03 H2-H2 
total CCSD(T) energies for these configurations are available from the EPAPS depository.— 
We have checked a small subset of these energies against calculations using the Dalton 
ab initio code;^^ the dimer total energies computed using the two codes agree to within 
2 X 10^® hartrees or better. 

To assess the error introduced by truncating the one-electron basis set at the aug-cc- 
pVQZ + (3s3p2d) level, we performed some calculations at selected configurations using a 
smaller aug-cc-pVTZ atom-centered basis set and the same (3s3p2d) bond function set. The 
coefficients Ai^i^^l obtained from these two sets of ab initio interaction energies are listed in 
Table HTl The two sets of coefficients generally differ by no more than 1% to 2%, suggesting 
that the aug-cc-pVQZ + (3s3p2d) basis set is nearly saturated. Truncating the one-electron 
basis set seems to have the largest effect on the isotropic coefficient Ao,o,o computed at small 
values of i?, where the potential energy surface is strongly repulsive. 

C. CCSDT ab initio calculations 

Our earlier study of the vibrational dependence of the H2-H2 interaction^ indicates 
that incompleteness in the many-electron basis set could materially affect the shape of the 
potential energy surface in the van der Waals well. Similar effects have been observed in other 
weakly-bound dimers of two-valence-electron systems.— 1^ To reduce the error associated 
with truncation of the many-electron basis set at the CCSD(T) level of theory, we employ 
a coupled- cluster treatment that includes a fully iterative treatment of single, double, and 
triple excitations,—"^ abbreviated CCSDT, to compute the H2-H2 interaction energy at 



selected high-symmetry geometries (those in which ^i, 6'2, and take values of or 7r/2). 
These calculations are performed using the tensor contraction engine^ incorporated into 
version 4.7 of the electronic structure code NWChem.— i^ 

Unfortunately, the CCSDT calculations are prohibitively expensive if we employ the 
aug-cc-pVQZ + (3s3p2d) one-electron basis set used in the CCSD(T) calculations. We 
therefore perform the CCSDT calculations using a smaller one-electron basis set consisting 
of only atom-centered aug-cc-pVTZ functions. We also use NWChem to perform CCSD(T) 
calculations at these high-symmetry geometries using the atom-centered aug-cc-pVTZ basis 
set. We then take the difference between the CCSDT and CCSD(T) counterpoise-corrected 
interaction energies as an additive correction to the aug-cc-pVQZ + (3s3p2d) CCSD(T) 
potential energy surface. For the sake of brevity, we will call this the "full-triples" correction. 
We found that to insure convergence of the CCSDT iterations at some geometries, it was 
necessary to increase the cutoff for computational linear dependence in the one-electron 
basis set to 10~^. For consistency, we therefore used this cutoff in all of the CCSDT and 
CCSD(T) calculations performed with NWChem. 

Because we compute the full-triples correction at a small number of H2-H2 orientations 
{61,62,4'), we cannot use the quadrature scheme described in the the previous subsection 
to extract corresponding full-triples corrections to the Ai^^i^^ coefficients computed at the 
CCSD(T) aug-cc-pVQZ + (3s3p2d) level of theory. Instead, we use least-squares techniques 
to fit the full-triples correction to the function 

AAo,o,o(i?, r-i, r2)Go,o,o(^i, 62, 4>) + AAo,2,2(i?, n, r2)G'o,2,2(^i, 62, 0) 

+ AA2,o,2(i?, ri, r2)G'2,o,2(^i, 62, <p) + AA2,2,4(^, ri, r2)G2,2,4(^i, 62, <p) . (6) 

We then add the corrections AA;^ jj ^ to the corresponding coefficients Ai^^i^l obtained 
from four-dimensional quadrature over the CCSD(T) aug-cc-pVQZ + (3s3p2d) interaction 
energies. The CCSDT and CCSD(T) energies used to compute the full-triples correction are 
available through EPAPS.'^^ For the sake of brevity, we henceforth use the term "coefficients" 
to mean the sum of the CCSD(T) coefficients and the full-triples corrections. 



D. Construction of a smooth potential energy surface 

We now construct a smooth potential energy surface from the ab initio coefficients 
Aij^,i2,L{R,ri,r2)- For each pair of H2 bond lengths (ri,r2), we create four cubic splines, 
one for each of the coefficients Ai-^^i^^l, that interpolate the 19 coefficient values between 
R = 4.25 Oo and R = 12.0 Oq. We extrapolate the splines to R values below 4.25 ag and 
above 12.0 Oq using functions described in the next two paragraphs. At R = 4.25 Qq, the 
slope of each cubic spline is constrained to match the slope of the corresponding small- i? 
extrapolating function. 

We extend each cubic spline to R values below 4.25 ao using a simple two-parameter 
exponential extrapolation of the form U exp{—cR) that fits the coefficients obtained at R = 
4.25 ao and 4.5 Oq. We should stress that this extrapolation is not expected to give highly 
accurate interaction energies for small R; we use it simply to define the slope for the cubic 
spline at i? = 4.25 Oq. The dimer bound state wave functions we compute using our potential 
energy surface are not sensitive to the highly repulsive small- i? region of the potential energy 
surface. 

Beyond i? = 12.0 ao, we extrapolate each spline using an inverse-power expansion of the 
form ^^ Cn/R^, including terms with n = 5 and 6 in the extrapolations for ^2,2.4, terms 
with 77, = 6, 8, and 10 for Aqq^q, and terms with n = 6 and 8 in the extrapolations for Aq^2,2 
and ^2,0,2- AH Cn coefficients are determined as functions of ri and ro. The C5 coefficient 



for ^2,2,4 is computed from the H2 quadrupole moments listed in Ref. 



47. The Ca coefficients 



are obtained from the isotropic and anisotropic R ^ dispersion energy coefficients given in 



Ref. 



48l and the expressions given in Ref. |49|. The Cg and Cio dispersion energy coefficients 



are obtained from Ref. 



50. 



To reduce the discontinuities in the higher-order derivatives of the coefficients at -R = 
12.0 ao, where the cubic spline meets the long-range inverse-power extrapolating function, we 
use the long-range function to compute values of the Ai^^i^^l coeffi cents at six evenly spaced 
"phantom" points ranging from R = 13.0 ao to -R = 18.0 ao, and force the spline to intercept 
these phantom points as well as the points computed at the 19 R values cited above. At 
R = 18.0 ao, we also constrain the slope of the spline to match that of the inverse-power 
expansion. However, we only use the spline to evaluate the coefficients between i? = 4.25 ao 
and R = 12.0 ao; beyond R = 12.0 ao, we use the inverse-power expansion to compute the 
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coefficients Ai^i^^^. 

Using these extrapolated cubic splines, we can compute the coefficients Ai^i^i{R,ri,r2) 

at any R for the discrete pairs of H2 bond lengths (ri,r2) at which we performed the ab 

initio calculations described above. As the last step in defining a smooth potential energy 

surface, we fit these interpolated (or extrapolated) coefficients to the expression 

2 2 
Y, H Cfc,„(ri - r,^)\r2 - r,^Y (7) 

fc=0 n=0 

where r^^ = 1.4 Oq. 

Figure [H shows how the isotropic coefficient AQfi^Q{R,ri,r2), vibrationally averaged over 
the ground state vibrational wave functions of the two H2 monomers, depends on R both 
in the small- i?, repulsive region of the potential energy surface and in the shallow H2- 
H2 well. We compare the vibrationally-averaged Aqqq coefficient computed in this work 
with a modified ab initio potential energy surface^ that gives accurate predictions for the 
low-temperature second virial coefficient of H2 gas, with the extrapolated CCSD(T) ab 
initio calculations of Diep and Johnson,—*^ and with an empirical isotropic potential energy 
curved obtained from an analysis of the total scattering cross section of moderate energy 
H2-D2 collisions. 

In the shallow well, our vibrationally averaged Ao,o,o coefficient agrees fairly well with 
the extrapolated CCSD(T) results,— >i^ which were computed within the rigid-rotor approx- 
imation using the v = vibrationally averaged bond length for both H2 monomers. The 
repulsive wall of our isotropic potential energy curve is slightly softer than that of the ex- 
trapolated CCSD(T) curve; our repulsive wall closely tracks the shape of the modified ab 
initio potential energy surface^^ that gives accurate second virial coefficients, except that 
our repulsive wall is shifted to slightly larger R values. It is interesting to note that in 
the small- i? repulsive region, the empirical isotropic potential energy curved derived from 
scattering data is considerably softer than any of the three curves derived from ab initio 
computations. 

III. COMPUTATION OF DIMER BOUND STATE ENERGIES 

We assess the quality of our potential energy surface by using it to compute the energies 
of several bound (and long-lived quasibound) states of the (112)2, H2-D2, and (02)2 dimers. 
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In this section, we summarize the methods used to compute these energies; the energies 
themselves are presented in later sections. We employ a standard close-coupled approach^- 
in which the nine-dimensional dimer wave function is written as 

^(R,ri,r2) = R-^J2^>^(^)^J'M,yi'R,h,r2)K,niri)(j)v,,nir2) ■ (§) 

A 

Here Tj is the vector separating the two nuclei of H2 molecule i, R is the vector separating 
the two molecules' centers of mass, and r^ and R are the corresponding unit vectors. The 
quantum numbers J and M represent respectively the total angular momentum of the dimer 
(excluding nuclear spin angular momentum) and its projection on a space-fixed z axis. The 
angular basis functions /j,m,7, which are defined as 

Ij,M,'yi'R,ri,r2) = ^ C{ji,mi,J2,fn2;Ji2,Mi2) 

X C(Ji2,Mi2,L,iV; J,M)F,,,^,(ri)r,,,™,(f2)ri,jv(R) , (9) 

couple the rotational angular momenta (ji, J2) of the two H2 molecules with the orbital 
angular momentum L of the dimer to create functions of definite J and M; we use 7 to 
represent the collection of angular momentum quantum numbers (ji, J2, Ji2i L), where J12 is 
the quantum number corresponding to the (vector) sum of the rotational angular momenta 
of the two H2 molecules. The summation index A represents a collection of eight quantum 
numbers: the four quantum numbers listed in 7, the total angular momentum quantum 
numbers J and M, and the vibrational quantum numbers Vi and f 2 of the two monomers. 
The functions (pvjif) are H2 monomer radial wave functions, defined so that 

Cj(^)0f',j'('") dr = 5^,v'5jj' (10) 

and obtained from a Numerov-Cooley^ analysis of the Kolos-Wolniewicz^ H2 potential 
energy curve. 

The dimer radial functions Fx{R) are solutions to a set of coupled second-order differential 
equations; the /^-dependent terms that couple the radial functions Fx{R) and Fxi{R) are 
obtained by integrating the six-dimensional H2-H2 potential energy surface over the eight 
coordinates (ri,r2,R), and are defined by replacing the rigid- rotor potential coefficients 



Ai^j2,LiR) in Eq- (9) of Ref. l30| with the corresponding vibrationally averaged coefficients 

/ / €ujM)K2,J2(^2)Ai,,i^,LiR,ri,r2)(l)v,,niri)(f)v2,J2ir2)dndr2 . (11) 

Jo Jo 
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We use the ABVN programa^ to evaluate the angular momentuin coupling coefficients that 



appear in Eq. (9) of Ref. 



3G. 



We convert the set of coupled second-order differential equations to a matrix eigenproblem 
by discretizing the equations on a grid in _R, ranging from -Rmin = 3.0 oq to -Rmax = 28.0 oq 
in steps of 0.1 ao, and replacing the dimer radial kinetic energy operator with a five-point 
central difference approximation evaluated on the grid. (Convergence tests show that using 
a smaller step size or a larger value of -Rmax does not change significantly the energies of the 
dimer states considered here.) We then solve the matrix eigenproblem using the ARPACK 
code^ driven by the SYMMLQ linear algebra routine.— We truncate the wave function 
expansion given in Eq. ([8]) by limiting ji and J2 to the values 0, 2, and 4, with the additional 
restriction that ji +J2 < 6, and by limiting vi and V2 to the values and 1. We also assume 
that the three vibrational manifolds defined by ft = fi -|- f2 = 0, 1, or 2 are effectively 
decoupled from one another, which further reduces the size of the matrix eigenproblem. 
The energies of the {v,j) rovibrational states of the H2 and D2 monomers and the dimer 
reduced masses, which appear in the close-coupled equations for the radial functions Fx{R), 
are computed from the parameters listed in Table IIIII 

Because we consider only even values of ji and J2 here, the parity of the angular basis 
function /j,a/,7 is controlled by the dimer orbital angular momentum quantum number L; 
when L is even, Ij^m,^ has even parity. Angular basis functions with different parities are not 
coupled together by Eq. ([8]). In addition, for a dimer of two identical monomers, the overall 
spatial wave function (exclusive of spin) must be either symmetric or antisymmetric under 
exchange of the two monomers, and the overall spin wave function must also be symmetric 
or antisymmetric under monomer exchange. The total wave function, which is the product 
of the spatial and spin wave functions, must be symmetric or antisymmetric under monomer 
exchange for bosonic and fermionic monomers, respectively. 

The para-B.2 molecule is a spin-zero composite boson. For a dimer of such bosons, no 
exchange-antisymmetric spin wave function can be constructed, and therefore only states 
whose spatial wave functions are symmetric under monomer exchange are physically admis- 
sible. These exchange-symmetric spatial wave functions are the only (112)2 wave functions 
considered here. On the other hand, ortho-T^2 molecules may have a total nuclear spin 
quantum number of either zero or two, and it is possible to construct {ortho-D2)2 dimers 
that have either an exchange-symmetric or an exchange-antisymmetric spin wave function. 
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Consequently the spatial wave function for {ortho-D2)2 may also be either symmetric or 
antisymmetric under monomer exchange, provided that the total (spin times spatial) (or- 
tho-D2)2 wave function is symmetric under monomer exchange.— 

To check that our matrix-based implementation of the close-coupled formalism is correct, 
we have used the BOUND code^^ to compute the energies of the (H2)2, H2-D2, and (02)2 
bound states that correlate with the monomers' {v,j) = (0,0) ground rovibrational states, 
and compare these energies with those obtained from our matrix-based code. [Because the 
BOUND code employs the rigid-rotor approximation, for this comparison we ignore the j 
dependence of the monomer radial wave functions 0t,,j(r) that appear in Eq. ([8]) and replace 
these radial wave functions with those for the monomers' ground rovibrational states. This is 
equivalent to neglecting centrifugal distortion effects on the monomer radial wave functions.] 
The good agreement between these two calculations confirms the validity of our matrix-based 
close-coupled approach. 

Some of the dimer states discussed below are long-lived quasibound states that can decay 
via rotational predissociation. The energies reported for these states are those obtained 
following the "infinite wall" procedure outlined by Grabenstetter and Le Roy,™ in which 
the energy of the quasibound state is monitored as -Rmax is decreased in 0.1 a^ steps. We 
estimate that using a finite step size of 0.1 oq in this procedure introduces an uncertainty in 
the quasibound state energies of no more than 0.002 cm~^. 

IV. EMPIRICAL ADJUSTMENTS TO THE POTENTIAL ENERGY SURFACE 

In this section, we show that if we make two small empirical modifications to our ab 
initio potential energy surface, it gives rotational and rovibrational transition energies for 
(112)2, H2-D2, and (02)2 dimers in good agreement with those obtained experimentally. The 
two modifications involve a small inward shift of the repulsive wall of the potential energy 
surface, which we quantify using an adjustable parameter s, and a slight increase in the 
magnitude of the surface's ^2,2,4 term, which we quantify using an adjustable parameter q. 
The unmodified, purely ab initio potential energy surface is defined by (s, q) = (0, 0). 

We focus first on H2-D2 and (02)2 dimer states that correlate with rotationally cold 
(j = 0) monomers as i? ^ 00. Because the wave functions of these states are overwhelmingly 
dominated by angular basis functions /j,m,7 with ji = j 2 = in Eq. (IHl), the states' energies 
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are insensitive to the anisotropic terms (^0,2,2, ^2,0,2, and ^2,2,4) of the potential energy 
surface; however, the states' energies are very sensitive to the location of the surface's 
repulsive wall. We therefore find the optimal value for s by adjusting s to bring the computed 
energies for transitions involving these states into good agreement with experimentally- 
measured transition energies. 

Next we consider IR-active transitions of the H2-D2 dimer which involve either (1) a 
pure vibrational transition f = 1 ^ in the H2 monomer and a pure rotational transition 
j = 2 ^ in the D2 monomer, or (2) a rovibrational transition (f , j) = (1, 2) <— (0, 0) in the 
H2 monomer and no excitation of the D2 monomer. These transitions involve final states 
whose energies are sensitive to the v4o,2,2 and ^2.0,2 anisotropic terms of the dimer potential 
energy surface. We find that, once the repulsive wall of the potential energy surface has been 
shifted inward slightly, the energies computed for these transitions are in good agreement 
with experimental measurements. This suggests that the ^0,2, 2 and ^2,0.2 terms of the shifted 
potential energy surface are accurate, at least in the range of R values probed by the H2-D2 
dimer wave functions. 

Finally, we examine (132)2 and (02)2 dimer states which correlate with R ^ 00 limits 
involving one j = and one j = 2 molecule. Some of these states have energies that are very 
sensitive to the strength of the ^2,2,4 term of the potential energy surface. By examining 
how the computed energies for transitions involving these states change with q, we find the 
value for q that gives the best overall agreement with experimental measurements. 

A. Combination differences from the H2 D2 and (02)2 dimer Qi(0) infrared spectra 

We begin by computing the J = 2 ^- spacings for the H2-D2 and (02)2 dimers that 
correlate with rotationally cold (j = 0) monomers; we perform these computations both for 
the dimers' f t = ground vibrational manifolds and for the f t = 1 manifold accessed by IR 
excitation of the H2 monomer in the H2-D2 dimer. Accurate experimental values for these 
J = 2 ^- spacings have been obtained from a combination-differences analysis of high- 
resolution H2-D2 and (02)2 IR absorption spectra.— Because of the large energy mismatch 
between the f = 1 levels of H2 and D2, in our calculations we assume that the H2-D2 dimer 
states correlating with H2 (f = 1) + D2 {v = 0) are decoupled from those correlating with 
H2 {v = 0) + D2 (f = 1). (Strictly speaking, the f t = 1 dimer states accessed in the IR 
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absorption experiment are quasibound, and can decay through vibrational predissociation. 
Our assumption that these states are decoupled from the f t = states, however, closes 
off this decay channel. Because the lifetimes of the f t = 1 dimer states are known to be 
extremely long,— this should not materially affect our results.) 

Table IIVI shows that the computed 7 = 2^—0 spacings are 0.015 cm~^ to 0.025 cm~^ 
lower than the experimental ones. If the dimers were rigid rotors, the J = 2 ^- spacings 
would be equal to six times the dimers' respective rotational constants. Because the dimers 
undergo large-amplitude zero-point motion along the R direction, a rigid-rotor model for 
the dimers' overall end-over-end rotation is not really appropriate. Nonetheless, this simple- 
minded picture suggests that the dimer states supported by the computed potential energy 
surface have average intermolecular distances that are slightly too large, by about 0.02 ao 
for the (02)2 dimer and 0.03 oq for the H2-D2 dimer. 

As we noted in our discussion of Fig. [H a small inward shift of the repulsive wall of 
our potential energy surface would bring it into closer agreement with a surface^ that 
gives accurate second virial coefficients for low-temperature H2 gas; such a shift would also 
reduce slightly the average intermolecular distances of the H2-D2 and (02)2 dimers, possibly 
bringing the computed J = 2 ^ spacings into better agreement with experiment. (We 



note here that the potential energy surface presented in Ref. Il8| was itself obtained by a 
similar inward shift of the repulsive wall of an ab initio computed potential energy surface.) 
We therefore modify our ab initio H2-H2 potential energy surface as follows. For R values 
below 6.5 ao, we shift our computed ab initio interaction energies to new, smaller, R values 
defined by 

-Knew = -Rold — ■S(6.5 flo " -Rold) , (12) 

and then construct a smooth s-dependent H2-H2 potential energy surface (as described 



above in Sec. I11D|) using the shifted points. Because we have not yet changed the strength 



of the ^2,2,4 term, we are at present implicitly holding q fixed at g = 0. 

Figure [2] shows how the errors in the J = 2 ^— spacings computed for the H2-D2 and 
(02)2 dimers change as s increases from s = to s = 0.025. Choosing s = 0.0175 brings 
all three of these computed spacings into in excellent agreement with experiment. Fixing s 
at this value amounts to an inward shift of the crossing point Rq, where the vibrationally- 
averaged H2-H2 isotropic coefficient Aq^qq{R) = 0, from Rq = 5.775 Qq for the original ab 
initio potential energy surface (with s = 0) to i?o = 5.762 ao for the empirically modified 
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surface. The corresponding shift for the H2-D2 dimer is from _Ro = 5.773 Oq to Rq = 5.760 Oq; 
for the (02)2 dimer, the shift is from Rq = 5.769 a^ to Rq = 5.756 a^. 

B. Qi(0) + So(0) and Si(0) infrared spectra of the H2 D2 dimer 

Next we consider transitions of the H2-D2 dimer in which either (1) the H2 monomer un- 
dergoes a pure f = 1 ^ vibrational transition and the D2 monomer simultaneously makes 
a pure j = 2 ^ rotational transition, or (2) the H2 monomer undergoes the rovibrational 
transition {v,j) = (1, 2) <— (0, 0) while the D2 monomer remains in its rovibrational ground 
state. The former transitions belong to the dimer's Qi(0) [H2] + So(0) [D2] band, and the 
latter transitions to the dimer's Si(0) [H2] band; for brevity, in this subsection we henceforth 
drop the molecular labels in square brackets and simply refer either Qi(0) + So(0) or Si(0) 
transitions. 

Because of the large energy mismatch between the j = 2 states of the H2 and D2 molecules, 
the upper states involved in these transitions are ones in which the j = 2 excitation remains 
localized on one of the monomers, and thus have energies that are insensitive to the ^2,2,4 
term of the potential energy surface. The computed energies for these transitions conse- 
quently provide insight into the quality of the surface's ^0,0,0, ^0,2,2, and ^2,0,2 terms. 

Four relatively sharp Qi(0) + So(0) transitions and two relatively sharp 5*1(0) transitions 
have been observed in the IR absorption spectrum of the H2-D2 dimer.— As Table |V] shows, 
the potential energy surface with s = 0.0175 gives transition energies for these six transitions 
in very good agreement with experiment; this suggests that the surface's v4o,2,2 and ^2,0,2 
terms are fairly accurate, at least over the range of R values for which the H2-D2 dimer has 
substantial probability density. 

C. So(0) infrared spectra of the (H2)2 and (02)2 dimers 

Finally we consider IR-active transitions of the (112)2 and (02)2 dimers that correlate 
with the So(0) j = 2 ^ pure rotational transitions of the H2 and D2 monomers. The 
upper states involved in these transitions are ones in which the j = 2 excitation is shared 
by the two monomers; the energies of these states are therefore sensitive to the ^2,2,4 term 
of the dimer potential energy surface, which couples together angular functions in Eq. ([8]) 
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with (ji,j2) = (0,2) and (ji,j2) = (2,0). 

Because of the low reduced mass of the (H2)2 dimer and the restrictions imposed by 
nuclear spin statistics, there is just one sharp So(0) IR-active transition for this dimer; it is 
a {J-iL) = (1,1) <— (0,0) transition and appears in the (H2)2 far-IR absorption spectrum 
at 355.425 cm~^.^- The transition energy computed for this absorption feature using the 
(s, g) = (0.0175,0) potential energy surface is 355.438 cm~^, or 0.013 cm"^ too high. 

The dimer wave function [Eq. ([8])] for the upper state of this transition contains signif- 
icant contributions from only four channels: those with {ji,J2,Ji2,L) angular momentum 
quantum numbers of (0, 2, 2, 1), (0, 2, 2, 3), (2, 0, 2, 1), and (2, 0, 2, 3). Furthermore, only two 
of these channels give independent contributions to the wave function; for the (H2)2 dimer, 
exchange symmetry constraints force the channels with {ji,J2,Ji2,L) = {a,b, Ji2,L) and 
(ji, J2, Ji2,L) = {b,a, Ji2,L) to, when L is odd, have radial functions Fx{R) that are equal 
in magnitude but opposite in sign. Figure [3] shows the Fx{R) radial functions for the two 
independent channels {ji,J2,Ji2,L) = (0,2,2,1) and (0,2,2,3) that define the upper state 
of the (J, L) = (1, 1) <— (0, 0) transition; about 97% of the upper state's probability density 
is associated with the two L = 1 channels. 

If we compute the expectation value of the dimer's potential energy using the upper state 
wave function, 

(^) = y y /"|^(R,ri,r2)|V(R,ri,r2)dRdridr2, (13) 

we find that it includes substantial contributions from the isotropic Ao,o,o term of the po- 
tential surface and the anisotropic ^0,2,2 and ^2,0,2 terms, along with a small contribution 
from the ^2^2,4 term; this last contribution is proportional to the integral 

/ FxiR)FxiR)dR (14) 

Jo 

where Fx and Fy are the two radial functions shown in Fig. [31 As we explained previously, 

the lower state for this transition has a wave function dominated by the (ji,J2) = (0,0) 

channel, and its potential energy expectation value is therefore sensitive to only the isotropic 

Ao,o,o term. 

This analysis suggests that a small perturbation of the ^2,2,4 term will change the energy 

of the upper state, but not that of the lower state, and could thus bring the computed 

transition energy for this far-IR absorption feature into better agreement with experiment. 

Furthermore, because the transitions considered in the preceding two subsections involve 
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states whose energies are insensitive to ^2,2,4, such a perturbation would preserve the good 
agreement with experiment observed for those transitions. (Naturally, we could also change 
the computed transition energy for this particular far-IR absorption feature by adjusting 
the Ao,2,2 and ^2,0,2 terms in the potential energy surface; however, such an adjustment 
would have the undesirable side effect of changing the transition energies computed in the 
immediately preceding subsection.) Here we adopt a very simple adjustment of the ^2,2,4 
term, which helps compensate for the fact that the ^2,2,4 coefficients computed in Sec. IIIAI 
include unwanted contributions from the electrostatic QH interaction: we multiply the ^2,2,4 
coefficients computed at each of the 19 R values by the quantity (1 + g), where q is an 
adjustable parameter, and then reconstruct the entire potential energy surface as described 
in Sec. HTPl 

Figure H] shows how q changes the computed position of the (H2)2 dimer's {J,L) = 
(1,1) <— (0,0) So(0) far-IR absorption feature. At g = 0.0235, the computed transition 
energy coincides with the experimental value of 355.425 cm^^. (However, the uncertainty of 
±0.005 cm~^ in this experimental transition energy means that a wide range of q values would 
be compatible with the experimental observations.) This suggests that a simple rescaling 
of our ^2,2,4 coefficients removes much of the QH interaction's erroneous contribution to 
these coefficients, even though the QH interaction has a different power-law dependence on 
R than does the quadrupole-quadrupole interaction that dominates the ^2,2,4 term. 

To place tighter constraints on q, we turn to the (02)2 dimer, which, because it is both 
heavier than (H2)2 and has less severe restrictions arising from nuclear spin statistics, exhibits 
many more absorption features in its far-IR So(0) band.— Twelve of these features are 
relatively sharp, suggesting that they involve bound or long-lived quasibound states, and 
also have firmly assigned initial- and final-state angular momentum quantum numbers. (We 
discuss later a thirteenth sharp transition whose initial- and final-state assignments are more 
tentative.) Figure [5] and Table |VT] show how the errors in the energies computed for these 
twelve transitions depend on q. For the g = potential energy surface, the deviations 
between computed and measured transition energies range from —0.025 cm^^ (transition e) 
to +0.064 cm~^ (transition h); at g = 0.0235, however, the computed energies for eight of the 
twelve transitions (those labeled e through /) agree with experiment to within ±0.007 cm^^. 
Only one of these eight transitions has a computed transition energy that differs from the 
experimental value by more than 0.005 cm~^, which is the experimental uncertainty quoted 
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for these transitions in Ref. 



28l . Furthermore, the value q = 0.0235 minimizes the mean 
absolute deviation between the predicted and observed transition energies for these eight 
transitions. 

The four transitions labeled a through d in Table |Vl] exhibit poorer agreement with 
experiment; furthermore, the errors in the computed energies for these four transitions at g = 
0.0235 are equal to or larger in magnitude than the errors at g = 0. Transition c corresponds 
to a very weak far-IR absorption feature; inspection of Fig. 3 in Ref. |28| shows that its 
intensity is comparable to the level of background noise in the (02)2 absorption spectrum, 
and it is possible that the true position of this feature differs slightly from that reported 



in Ref. 



28 



Transitions a, b, and d, however, correspond to relatively strong absorption 
features; furthermore, while transition a is a shoulder on the low-energy side of a very intense 
feature (transition /), transitions b and d are well isolated from other spectral features, and 
transitions a and b are linked by the J = 3 ^ 1 combination difference of the dimer's 
Ji = J2 = manifold. It seems unlikely that the quoted experimental uncertainties for these 
three transitions could be badly underestimated. It therefore appears that our potential 
energy surface slightly underpredicts the energies of the (J, L) = (2, 2) and (3, 3) excited- 
state levels accessed via these three IR transitions. 

The experimental (02)2 So(0) IR absorption spectrum exhibits a thirteenth sharp feature, 
corresponding to the transition energy 176.627 cm~^, which might be either the {J,L) = 
(3, 1) <— (2, 2) transition or the (J, L) = (3, 2) ^- (3, 3) transition.— Our potential energy 
surface predicts transition energies of 176.645 cm~^ and 176.598 cm~^, respectively, for these 
transitions. 

Finally, we note that the computed transition energies listed in Table |V] change by only 
0.001 to 0.002 cm~^ when the {s, q) = (0.0175, 0.0235) potential energy surface is used. This 
validates our decision to hold q fixed at g = while we find the optimal value for s, and 
then hold s fixed at this value while we find the optimal value for q. 

V. OTHER COMPARISONS WITH EXPERIMENT 

In the previous section, we showed that the quality of the four Ai-^j^.l terms of our 
potential energy surface could be assessed individually by considering transitions between 
pairs of states that have energies sensitive to specific subsets of these terms. We found that 
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with two small adjustments to the potential energy surface, we could generate a surface that 
gives computed transition energies in fairly good agreement with a number of high-resolution 
experimental measurements. 

Although some of the transitions considered in the previous section involve vibrational 
excitation of the H2 monomer in the H2-D2 dimer, we have not yet considered vibrationally 
excited states of the (H2)2 or (02)2 dimers. In these dimers' Wt = 1 vibrationally excited 
states, the vibrational excitation is delocalized across the pair of monomers; transitions to 
these excited states therefore probe the simultaneous dependence of the potential energy 
surface on ri and r2. 

In this section, we show that our modified potential energy surface predicts energies for 
these transitions that are in good agreement with experiment, indicating that the surface 
accurately describes the vibrational coupling between the two monomers in the (112)2 and 
(02)2 dimers. We also consider IR-active double vibrational transitions of the (02)2 dimer, 
in which each monomer undergoes a f = 1 <— excitation; the good agreement we obtain 
with experiment provides further evidence that our modified potential energy surface has 
the correct (ri,r2) dependence. 

A. Qi(0) spectra of the (112)25 H2 D2, and (02)2 dimers 

Tables NTH through HXl list the energies of several (112)2, H2-D2, and (02)2 bound states 
that correlate with j = monomer states, computed using the final (s, q) = (0.0175, 0.0235) 
potential energy surface. Using these bound state energies and the monomer Qi(0) transition 
energies from Table IIIH we can obtain theoretical positions for the P and R lines in the 
dimers' Qi(0) IR absorption spectra. In Table 1X1 we list the computed positions for the 
eleven P and R lines that have been observed experimentally,'^ and compare the computed 
positions with the observed ones. 

The computed transition energies for (112)2 and H2-D2 are in excellent agreement with 
experiment, deviating from the observed energies by amounts smaller than the estimated ex- 
perimental uncertainties. The transition energies for (02)2, however, deviate systematically 
from the experimental measurements by about —0.01 cm^^, or about twice the estimated 
uncertainty in the measured transition energies. 

To investigate this discrepancy further, we have computed the transition energies of the 
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(02)2 dimer's P(2), P(l), R(0), and R(l) lines using a set of potential energy surfaces with 
different s values, keeping q fixed at g = 0.0235. (The P and R lines involving J = 3 states 
have been omitted from this analysis simply because computing these states' energies at 
several values of s is very time consuming.) In Fig. [6] we show how the deviations between 
the computed and observed transition energies change with s. Only for s values near 0.0175 
do the computed transition energies deviate systematically from experiment; in addition, 
the s = 0.0175 energies listed in Table [IX] give J = 2 ^ and 3^1 spacings for both 
the f t = and IR-active f t = 1 manifolds within 0.002 cm~^ of the experiment ally- derived 
values.— These observations suggest that the systematic deviations observed for (02)2 in 
Table |X] are not related to a poor choice for s. 

These discrepancies could indicate a small error in the isotropic Ao,o,o term's simultaneous 
dependence on ri and r2; the vibrationally excited (02)2 states involved in the transitions 
listed in Table 1X1 are antisymmetric linear combinations of (^1,^2) = (0, 1) and (1,0) states 
and are therefore sensitive to this (ri,r2) coupling. The same coupling term, however, 
is also active in the IR-active vibrationally excited (112)2 state, and in the Raman-active 
(112)2 excited state discussed in the next paragraph, and the agreement with experiment is 
excellent for transitions involving these states of the (112)2 dimer. More work is needed to 
understand the systematic deviations in the final column of Table IIX[ 

The (112)2, H2-D2, and (02)2 dimers should all have Raman-active transitions in the 
vicinity of the monomers' Qi(0) Raman transitions; thus far, however, only the (112)2 dimer's 
Raman spectrum has been observed experimentally.— It consists of a single narrow line 
corresponding to a transition energy of 4160.78 ± 0.02 cm^^. The theoretical Qi(0) Raman 
transition energy for (112)2 derived from the first two lines of Table IVnl is 4160.764 cm~^; 
the difference between the computed and experimental Raman transition energies is only 
slightly smaller than the estimated experimental uncertainty. However, the reported value 
for the experimental Raman transition energy is based on a value of 4161.18 cm~^ for the 
Qi(0) transition of an isolated H2 molecule. The difference between this value and the value 
used here (4161.169 cm~^) accounts for more than half the difference between the computed 
and observed dimer Raman transition energies. If we instead compare the observed and 
computed dimerization-induced red shift of the II2 Qi(0) Raman transition, we find that 
our computed red shift of 0.405 cm~^ is in excellent agreement with the reported value^^ of 
0.400 ± 0.02 cm-^ 
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B. Si(0) and Qi(0) + Si(0) infrared spectra of the (H2)2 and (02)2 dimers 

We finally use our modified potential energy surface to compute transition energies for 
features in the Si(0) IR absorption bands of the (H2)2 and (02)2 dimers, and for features 
in the Qi(0) + Si(0) IR absorption band of the (02)2 dimer. The (02)2 transitions involve 
upper states whose energies are sensitive to the (ri,r2) dependence of the potential energy 
surface. 

The Si(0) IR absorption spectrum of (112)2 contains just one narrow line,— at 4498.734 ± 
0.004 cm^^. This feature is associated with a transition from the dimer's ground state 
(the first line of Table |Vl) to a state with J = L = 1 that is a linear combination of 
(fi, ji; f2, J2) = (1,2; 0,0) and (0,0; 1,2); the final state's energy is listed in the first line 
of Table IXll The {s,q) = (0.0175,0.0235) potential energy surface gives a computed Si(0) 
dimer transition energy of 4498.729 cm~^, in excellent agreement with the observed value. 

The angular basis functions /j,m,7 in Eq- & that correspond to (J, L, ji, J2) = (1, 1, 0, 2) 
and (1, 1, 2, 0), which dominate the final state wave function for this Si(0) dimer transition, 
are not directly coupled together by any of the four terms A/i,/2,l that appear in our potential 
energy surface; consequently, this (112)2 transition, like the H2-D2 transitions considered 
in Sec. lIVBt probes primarily the monomer vibrational dependence of the surface's Aqq^q, 
^0,2,2; and v42,o,2 terms. In contrast to the H2-D2 transitions discussed in Sec. IIVB| however, 
the (112)2 Si(0) transition is sensitive to the ^2,2,0 and ^2,2,2 terms of the potential energy 
surface, which we have ignored; the fact that we obtain good agreement with experiment 
without explicitly including these terms in our surface is further evidence that these terms 
are of minor importance for the dimer bound states considered in this work. 

The Si 



in Ref. 
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_ 0) IR absorption band for (02)2 is much richer than that for (112)2, and is described 
as "possibly [the] most informative of all the hydrogen dimer spectra" presented 
there. It contains six pairs of narrow lines separated by the (02)2 f t = ground-state 
J = 2<— 0or3^-l spacings (3.001 cm~^ and 4.814 cm~^, respectively) and nine additional 
narrow lines. 

The six pairs of lines are associated with transitions from two different rotational levels of 
the f t = ground state to a common f t = 1 upper state level with a firm angular momentum 
quantum number assignment. The upper part of Table Kill compares the computed and 
experimental transition energies for the higher-frequency transition of each of these pairs. 
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(No additional information about the quality of our potential energy surface is carried by the 
other transition of each pair.) The agreement between computed and measured transition 
energies is quite satisfactory; the largest deviation is 0.030 cm~^, for the {J,L) = (2,4) <— 
(1, 1) transition. 
In Ref. 26 



initial- and final-state quantum number labels were proposed for the nine other 
narrow lines that appear in the (02)2 dimer's Si(0) IR absorption band; these assignments 
were descrbed as "less certain" than the assignments for the pairs of lines linked by ground- 
state combination differences. Eight of these lines are listed in the lower part of Table Kill 
which shows that using these transition assignments, we again observe very good agreement 
between computed and measured transition energies. (The one line omitted from Table IXIII 
involves a transition to a J = 4 state whose energy we have not attempted to compute.) 
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Table IXIII thus confirms the transition assignments proposed in Ref. 

The upper-state wave functions for the transitions listed in Table IXlI] are linear combina- 
tions of radial functions in Eq. ^ with angular momentum quantum numbers (J, L, ji, J2) = 
(J, L, 0,2) and (J, L, 2,0). For L > 2, these pairs of radial functions are coupled together 
by the ^2,2,4 term of the potential energy surface. In addition, some of the upper states 
accessed via these transitions are mixtures of Si(0) states, in which the rotational and vi- 
brational excitation reside on the same D2 monomer, and Qi(0) + So(0) states, in which 
one monomer is vibrationally excited while the other is rotationally excited; the quantity F 
listed in Table IXIII measures the degree of mixing in the upper-state wave functions. The 
good agreement between computed and observed transition energies in Table IXIII especially 
for transitions to upper states with L > 2 or with F values below 0.9, indicates that the 
(ri,r2) dependence of our potential energy surface, and of the ^2,2,4 term in particular, is 
reasonably accurate. 

Further evidence that the (ri,r2) dependence of our potential energy surface is accurate 
comes from Table IXIIIl where we compare the computed and observed transition energies 
for several (02)2 transitions in the dimer's Qi(0) + Si(0) overtone IR absorption band. The 
agreement between computed and observed transition energies is fairly good, although it 
appears that the potential energy surface generally underestimates slightly the energies of 
the upper states of these transitions. 
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VI. PREDICTIONS FOR NOT- YET-OBSERVED TRANSITIONS 

In this section, we use our final (s, q) = (0.0175, 0.0235) potential energy surface to 
predict the energies of some not-yet-observed IR transitions of the H2-D2 dimer and some 
not-yet-observed Raman transitions of the (112)2 and (02)2 dimers. Experimental studies 
designed to search for these transitions would help test the accuracy of the potential energy 
surface presented here. 

The region of the H2-D2 dimer's IR absorption spectrum associated with Qi(0) excita- 
tion of the H2 monomer has already been studied experimentally, and as Table 1X1 shows, our 
potential energy surface gives accurate transition energies for the four P and R lines in this 
portion of the dimer's IR spectrum. The dimer should have four additional IR-active P and 
R lines associated with Qi(0) excitation of the D2 monomer. We have computed the tran- 
sition energies for these four lines based on the dimer binding energies listed in Table IVIIH 
Table IXIVI lists the predicted transition energies for these four absorption features. 

As we noted earlier, the single Raman-active transition in the (112)2 dimer's Qi(0) 
band was recently observed. This dimer should also have Raman-active transitions in the 
monomer So(0) and Si(0) bands. The energies for these two transitions can be computed 
from the dimer's ground state binding energy of 2.894 cm~^ and the binding energies of 
the (J, L) = (2, 0) states listed in the second line of Table IXII We therefore predict that 
the (112)2 dimer should exhibit So(0) and Si(0) Raman transitions at 354.245 cm~^ and 
4497.415 cm~^, respectively. Unfortunately, these transitions are fairly close to the corre- 
sponding Raman-active transitions of the free H2 monomer, which are located at transition 
energies of 354.373 cm~^ and 4497.839 cm~^, so a high-resolution experiment will likely be 
required to observe the dimer transitions. 

Because the (02)2 dimer has four bound states, its Raman spectrum will be much richer 
than that of (112)2. In Tables IXVl and IXVII we give predictions for Raman- active transitions 
of the (02)2 dimer in the monomer Qi(0), So(0), and Si(0) bands; these predictions are 
based on the energy levels listed in Tables IIXI and IXVIII There will be additional Raman 
features in the dimer's So(0) and Si(0) bands, associated with transitions to final dimer 
states with J = 3, which we have omitted from Table IXVII because we have not computed 
the energies of these final dimer states. 
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VII. SUMMARY AND DISCUSSION 

We have presented a six-dimensional H2-H2 potential energy surface that accurately 
describes several bound (and long-lived quasibound) states of the (H2)2, (02)2, and H2-D2 
dimers that correlate with H2 and D2 monomers in their (f , j) = (0, 0), (0, 2), (1, 0), and (1, 2) 
rovibrational states. The surface is based on a set of ab initio H2-H2 interaction energies 
that appear to be nearly converged with respect to the one-electron and many-electron basis 
sets, and which cover fairly densely the region of configuration space associated with the 
dimer's van der Waals well. The surface incorporates two empirical adjustments: one softens 
slightly the surface's repulsive wall at small H2-H2 distances, and one increases slightly the 
magnitude of the surface's ^2,2,4 term that couples the rotational degrees of freedom of 
the two monomers. The latter adjustment appears to compensate for the fact that our 
original ^2,2,4 coefficients include small contributions from the electrostatic quadrupole- 
hexadecapole interaction between the two molecules. A Fortran subroutine that evaluates 
the {s, q) = (0.0175, 0.0235) potential energy surface is available through EPAPS.— 

An empirical softening of the ab initio H2-H2 interaction energy at small intermolecular 
distances is not unprecedented. For instance, Schaefer and Kohler— found that by softening 
the short-range repulsive wall of an ab initio H2-H2 potential energy surface, they could 
bring properties computed using the surface (in this case, the second virial coefficient of 
H2 gas) into better agreement with experiment. The surface of Diep and Johnson^^ii^ is 
based on a series of ab initio calculations that have been extrapolated to the estimated 
complete one-electron basis set limit; this extrapolation technique yields a surface whose 
repulsive wall is slightly softer than that of the largest basis-set surface explicitly calculated, 
and thus has the same effect as an empirical softening of the repulsive wall. It would be 
interesting to compute the H2-H2 interaction energy using explicitly-correlated electronic 
structure methods^ to see whether lingering basis set incompleteness in the present ab 
initio calculations is what necessitates the softening of the short-range repulsive wall of the 
surface. 

We have used our potential energy surface to predict the energies of 34 not-yet-observed 
IR and Raman transitions for (1^2)2, (02)2, and H2-D2 dimers involving even-j states of the 
H2 and D2 monomers. Observations of these transitions could help verify the accuracy of 
the present potential energy surface, or point out areas where the surface could be further 
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improved. Calculations of the energy levels of dimers containing one j = 1 molecule, such 
as ortho-}l2-para-B.2, could also be useful in this regard; these calculations are in progress 
and will be reported in due course. 
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FIG. 1: (Color online.) Comparison of the vibrationally-averaged isotropic H2-H2 potential energy 
curvce ^0,0,0 (-R) obtained in this work (solid line) with the isotropic rigid-rotor potential energy 
curves obtained by other researchers. Boxes represent the extrapolated CCSD(T) potential energy 



surface of Refs. 



14 and 



181 : the dashed 



15l : circles represent the adjusted ab initio potential of Ref. 
line represents the empirical potential of Ref. |5l|. The solid line shown here is computed from the 
unmodified (s, q) = (0, 0) potential energy surface. Panel (a) shows the repulsive wall at small R 



values; panel (b) shows the van der Waals well. 
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FIG. 2: (Color online.) Dependence on s of the errors (computed minus experiment) in the 
J = 2 ^ spacings of the H2-D2 and (02)2 dimers. Solid lines are for the vibrationally cold 
(ft = 0) dimers; the dashed line is for the H2 {v = 1) + D2 (v = 0) dimer. The parameter q is held 
fixed at g = 0. 
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FIG. 3: (Color online). The two dominant Fx(R) radial functions for the {J,L) = (1,1) upper 
state accessed in the (H2)2 dimer's So(0) IR transition. These functions are computed using the 
(s, q) = (0.0175, 0) potential energy surface. 
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FIG. 4: Dependence on q of the error (computed minus experiment) in the transition energy 
computed for the (H2)2 dimer's far-IR So(0) absorption feature. The parameter s is held fixed at 
s = 0.0175. 
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FIG. 5: (Color online.) Dependence on q of the errors (computed minus experiment) in the 
transition energies computed for twelve features in the (D2)2 dimer's far-IR So(0) absorption band. 
Panel (b) is a magnification of the small box in panel (a). The labels affixed to each line refer to 
Table IVTl The parameter s is held fixed at s = 0.0175. 
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FIG. 6: Dependence on s of the errors (computed minus experiment) in the transition energies 
computed for four hues in the (02)2 dimer's Qi(0) IR absorption band. The parameter q is held 
fixed atq = 0.0235. 
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TABLE I: Comparison of angular expansion coefficients Ai^^i^^i^ (in cm^^) computed using two 
spherical quadrature rules at (i?, ri, r2) = (4.5 ao, 1.4 ao) 1-7 ao). These coefficients do not include 
the full-triples correction. 

{li,l2,L) 18-point 24-point 

(0, 0, 0) 674.559 674.629 

(0, 2, 2) 32.984 33.577 

(2, 0, 2) 20.195 20.465 

(2, 2, 4) 19.017 19.174 

(2, 2, 0) -0.599 1.413 

(2, 2, 2) 0.465 -0.292 
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TABLE II: Angular expansion coefficients Ai^^i^^i^^R, ri,r2), in cm^-*^, computed from aug-cc-pVTZ 
and aug-cc-pVQZ ab initio energies. A (3s3p2d) set of bond functions is used in all calculations. 
The coefficients are grouped into pairs of rows corresponding to fixed (R, ri, r2); the upper row in 
each pair lists the aug-cc-pVQZ coefficients, while the lower row in each pair lists the aug-cc-pVTZ 
coefficients. These coefficients do not include the full-triples correction. 



Riao) 



iri,r2) (ao) 



iO,0,0 



A 



2,0,2 



A, 



0,2,2 



A 



2,2,4 



4.5 



4.5 



4.5 



4.5 



(1.1, 1.1) 


378.145 


10.026 


10.026 


6.360 




382.457 


10.094 


10.094 


6.320 


(1.4, 1.4) 


570.418 


19.989 


19.989 


14.277 




576.072 


20.116 


20.116 


14.197 


(1.7, 1.7) 


776.340 


31.141 


31.141 


24.607 




783.772 


31.379 


31.379 


24.478 


(1.1, 1.7) 


573.576 


11.959 


33.199 


13.032 




579.204 


11.878 


33.480 


13.044 


(1.1, 1.1) 


101.523 


3.023 


3.023 


3.607 




103.343 


3.046 


3.046 


3.585 


(1.4, 1.4) 


174.837 


6.794 


6.794 


8.119 




177.554 


6.879 


6.879 


8.070 


(1.7, 1.7) 


262.487 


11.883 


11.883 


14.347 




266.336 


12.043 


12.043 


14.254 


(1.1, 1.7) 


178.370 


3.853 


12.040 


7.370 




181.142 


3.883 


12.179 


7.325 



5.0 



5.0 



5.0 



5.0 



6.5 



(1.4, 1.4) 



-22.530 
-22.239 



-0.543 
-0.528 



-0.543 
-0.528 



2.070 
2.054 
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TABLE III: Monomer spectroscopic constants (in cm ^) and total masses employed in the dimer 
bound state calculations. 



Species 



AE (v = l 



0) 



B {v = 0) 



B (v = l) 



Mass (m-e) 



H2 
D2 



4161.169 
2993.614 



59.0622 

29.8445 



56.1117 

28.7908 



3674.3 
7342.9 



TABLE IV: Observed and computed spacings (in cm~^) between the J = and J = 2 states of 
H2-D2 and D2-D2 dimers that correlate with j = states of the constituent monomers. These 
bound state computations employ the unmodified (s, q) = (0, 0) potential energy surface. 



Dimer 



Observed 



Computed 



H2 (?; = 0) + D2 {v = 0) 
H2 (^ = 1) + D2 {v = 0) 
D2 {v = 0) + D2 {v = 0) 



3.848 
3.889 
3.001 



3.822 
3.866 
2.986 



TABLE V: Observed transition energies and deviations between observed and computed transition 
energies (computed minus experiment), in cm~^, for selected IR-active transitions of the II2-D2 
dimer involving v = 1 <— excitation of the II2 monomer. The D2 monomer remains in its w = 
vibrational level during the transition. The initial and final states are identified by the angular 
momentum quantum numbers associated with the dominant term in the dimer wave function 
[Eq. ([8])]; ji and j2 are the angular momenta of the II2 and D2 molecules, respectively. The 
computed transition energies are obtained using the (s,g) = (0.0175,0) potential energy surface. 



{J,L,n,J2y^iJ,L,n,J2)" 



Observed 



Deviation 



(1, 1, 0, 2) 
(2, 1, 0, 2) 
(0, 2, 0, 2) 
(1, 2, 0, 2) 
(1, 1, 2, 0) 
(2, 1, 2, 0) 



(2, 2, 0, 0) 

(2, 2, 0, 0) 

(1, 1, 0, 0) 

(1, 1, 0, 0) 

(2, 2, 0, 0) 

(2, 2, 0, 0) 



4337.046 
4337.609 
4342.004 
4342.208 

4494.719 
4495.20 



+0.003 
+0.002 
+0.013 
+0.007 
-0.007 
+0.009 
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TABLE VI: Observed transition energies and deviations between observed and computed transition 
energies (computed minus experiment), in cm^-*^, for IR-active (02)2 transitions in the D2 So(0) 
band. The computed transition energies are obtained using either the {s,q) = (0.0175,0) or the 
(■5) q) = (0.0175, 0.0235) potential energy surface. 



Label 



(J, ly ^ (J, L)" 



Observed 



Deviation 
q = 



Deviation 
q = 0.0235 



a 


(2, 


2). 


-(3 


3) 


175.507 


+0.010 


-0.012 


h 


(2, 


2). 


-(1 


1) 


180.322 


+0.005 


-0.016 


c 


(2, 


1)- 


-(2 


2) 


177.359 


-0.022 


-0.026 


d 


(3, 


3)- 


-(2 


2) 


181.287 


-0.010 


-0.029 


e 


(1, 


2). 


-(1 


1) 


182.328 


-0.025 


+0.007 


/ 


(1, 


3)- 


-(0 


0) 


184.536 


-0.006 


+0.005 


9 


(2, 


0)- 


-(3 


3) 


172.776 


+0.024 


+0.004 


h 


(0, 


2). 


-(1 


1) 


177.996 


+0.064 


+0.000 


i 


(2, 


3)- 


-(2 


2) 


182.797 


-0.024 


-0.002 


J 


(2, 


0)- 


-(1 


1) 


177.592 


+0.018 


-0.001 


k 


(1, 


1)- 


-(0 


0) 


178.747 


+0.043 


+0.002 


I 


(1, 


1)- 


-(2 


2) 


175.744 


+0.045 


+0.004 



TABLE VII: Binding energies (in cm^"*^) of selected (112)2 dimer states that correlate with 
J = states of the constituent monomers. These bound state computations employ the 
(S)?) = (0.0175,0.0235) potential energy surface. 



Vt 



J 



Binding energy 




1 
1 





1 



2.894 
3.299 
1.554 
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TABLE VIII: Binding energies (in cm^^) of selected H2-D2 dinier states that correlate with 
j = states of the constituent monomers. These bound state computations employ the 
(S)?) = (0.0175,0.0235) potential energy surface. 

J Vibrational state 

?;t = v = l (H2) V = 1 (D2) 

4.417 4.792 4.644 

1 3.074 3.442 3.299 

2 0.568 0.902 0.775 

TABLE IX: Binding energies (in cm~^) of selected (02)2 dimer states that correlate with j = 
states of the constituent monomers. The letters S and A indicate that the state is respectively 
symmetric or antisymmetric under monomer exchange. These bound state computations employ 
the {s, q) = (0.0175, 0.0235) potential energy surface. 

J Vibrational state 

«t = «t = 1 (S) i;t = 1 (A) 

6.712 (S) 7.118 6.939 

1 5.696 (A) 5.925 6.101 

2 3.711 (S) 4.111 3.941 

3 0.885 (A) 1.109 1.263 



40 



TABLE X: Computed transition energies (in cm^^) and deviations from experiment (computed mi- 
nus experiment) of the P and R lines in the Qi(0) IR bands of the (H2)2, H2-D2, and (02)2 dimers. 
For the H2-D2 dimer, the vibrationally excited state correlates with v = 1 II2 + v = D2; for 
the (112)2 and (02)2 dimers, the vibrational excitation is delocalized antisymmetrically across the 
two monomers. The computed transition energies are obtained using the (s, q) = (0.0175, 0.0235) 
potential energy surface. 



J'^J" (H2)2 H2-D2 (D: 



2i2 







Computed 


Deviation 


Computed 


Deviation 


Computed 


Deviation 


2 <- 


-3 


— 


— 


— 


— 


2990.558 


-0.007 


1 ^ 


-2 


— 


— 


4158.296 


-0.003 


2991.400 


-0.012 


0^ 


- 1 


— 


— 


4159.451 


-0.001 


2992.371 


-0.011 


1 ^ 


-0 


4162.509 


-0.004 


4162.144 


-0.003 


2994.401 


-0.011 


2^ 


-1 


— 


— 


4163.341 


+0.000 


2995.368 


-0.011 


3^ 


-2 


— 


— 


— 


— 


2996.216 


-0.013 



TABLE XI: Energies (in cm^^) for (112)2 states with J12 = 2, wt = or 1, and L < 1. The energies 

are obtained using the (s, q) = (0.0175, 0.0235) potential energy surface, and are given relative to 
the So(0) and Si(0) H2 monomer energies for f t = and 1, respectively. 

(J, L) Energy Energy 

■ut = wt = 1 

(1, 1) -1.842 -2.004 

(2, 0) -3.022 -3.318 

(2, 1) -1.191 -1.463 
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TABLE XII: Observed transition energies and deviations between observed and computed tran- 
sition energies (computed minus experiment), in cm^^, for IR-active (02)2 transitions in the D2 
Si(0) band. The computed transition energies are obtained using the (s,g) = (0.0175,0.0235) 
potential energy surface. The column labeled F indicates the fraction of upper-state probability 
associated with functions in Eq. ([8]) in which the v = 1 and j = 2 molecular excitations reside on 
the same monomer. 



(J, 


L)'^ 


-{J, 


L)" 


Observed 


Deviation 


F 


(2 


0)- 


-(1, 


1) 


3164.705 


-0.005 


0.984 


(1 


1)- 


-(0, 


0) 


3166.195 


-0.005 


0.959 


(2 


2). 


-(1, 


1) 


3167.890 


-0.004 


0.966 


(1 


3)- 


-(0, 


0) 


3169.919 


-0.020 


0.666 


(2 


0)- 


-(1, 


1) 


3170.429 


-0.020 


0.196 


(2 


4)- 


-(1, 


1) 


3173.702 


+0.030 


0.467 


(3 


2)- 


-(3, 


3) 


3163.378 


-0.009 


0.993 


(3 


1)- 


-(2, 


2) 


3163.707 


-0.003 


0.989 


(2 


1)- 


-(2, 


2) 


3164.281 


-0.011 


0.992 


(0 


2)- 


-(1, 


1) 


3166.340 


-0.004 


0.871 


(1 


2)- 


-(1, 


1) 


3167.392 


-0.030 


0.942 


(2 


3)- 


-(2, 


2) 


3168.343 


-0.028 


0.860 


(3 


3)- 


-(2, 


2) 


3168.737 


-0.004 


0.891 


(2 


1)- 


-(2, 


2) 


3170.931 


+0.003 


0.298 
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TABLE XIII: Observed transition energies and deviations between observed and computed transi- 
tion energies (computed minus experiment), in cm^^, for IR-active (02)2 transitions in the D2 Qi(0) 
+ Si(0) band. The computed transition energies are obtained using the {s,q) = (0.0175,0.0235) 
potential energy surface. 

(J, L)' ^ (J, L)" Observed Deviation 

6152.870 -0.018 

6155.672 -0.020 

6155.672 -0.021 

6156.767 -0.015 

6157.770 -0.003 

6158.669 -0.016 

6161.471 -0.036 

6162.773 -0.002 

6163.261 -0.015 

6164.872 -0.005 

TABLE XIV: Predicted transition energies (in cm^^) of the P and R hues in the Qi(0) IR band 
of the H2-D2 dimer, for vibrationally excited states correlating with t; = 0H2 + w = lD2. The 
predictions are obtained using the {s, q) = (0.0175, 0.0235) potential energy surface. 

J' ^ J" Energy 

1^2 2990.884 

^ 1 2992.044 

1 ^ 2994.732 

2 ^ 1 2995.913 



(2, 0) ^ 


- (3, 3) 


(1, 1) <- 


- (2, 2) 


(2, 2) . 


- (3, 3) 


(3, 1) . 


- (2, 2) 


(0, 2) . 


- (1, 1) 


(1, 1) <- 


- (0, 0) 


(3, 3) ^ 


- (2, 2) 


(1, 2) <- 


- (1, 1) 


(2, 3) ^ 


- (2, 2) 


(1, 3) . 


- (0, 0) 



(1, 1) - 


- (3, 3) 


(0, 0) <- 


- (2, 2) 


(0, 0) <- 


- (0, 0) 


(1, 1) - 


- (1, 1) 


(2, 2) <- 


- (2, 2) 


(3, 3) <- 


- (3, 3) 


(2, 2) . 


- (0, 0) 


(3, 3) ^ 


- (1, 1) 
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TABLE XV: Predicted transition energies (in cm^^) for Raman-active transitions in the Qi(0) 
band of the (D2)2 dimer. The predictions are obtained using the {s,q) = (0.0175,0.0235) potential 
energy surface. 

(J, ly ^ (J, Lf Energy 

2988.398 
2990.207 
2993.208 
2993.209 
2993.214 
2993.236 
2996.215 
2998.047 



TABLE XVI: Predicted transition energies (in cm~^) for Raman-active transitions in the So(0) and 
Si(0) bands of the (02)2 dimer. The predictions are obtained using the {s,q) = (0.0175,0.0235) 
potential energy surface. 

(J, LY ^ (J, L)" So(0) energy Si(0) energy 

183.837 3167.595 

178.713 3165.725 

182.694 3168.981 

178.108 3165.208 

184.568 3168.916 

180.836 3164.594 

175.712 3162.724 

179.693 3165.980 

173.298 3160.397 

179.758 3164.105 



(0, 2) ^ 


- (0, 0) 


(2, 0) . 


- (0, 0) 


(2, 2) <- 


- (0, 0) 


(1, 1) - 


- (1, 1) 


(1, 3) <- 


- (1, 1) 


(0, 2) . 


- (2, 2) 


(2, 0) . 


- (2, 2) 


(2, 2) . 


- (2, 2) 


(1, 1) - 


- (3, 3) 


(1, 3) <- 


- (3, 3) 



44 



TABLE XVII: Energies (in cm-^) for (02)2 states with J12 = 2, wt = or 1, and J < 2. The 
energies are obtained using the {s, q) = (0.0175, 0.0235) potential energy surface, and are given 
relative to the So(0) and Si(0) D2 monomer energies for f t = and 1, respectively. 

(J, L) Exchange Energy Energy 

symmetry ^'t = ^'t = 1 



(0 


2) 


S 


-1.942 


-5.476 


(1 


1) 


S 


-7.030 


-6.881 


(1 


2) 


s 


-5.620 


-4.850 


(1 


3) 


s 


-1.238 


-3.172 


(2 


0) 


s 


-7.066 


-7.347 


(2 


1) 


s 


-5.445 


-5.800 


(2 


2) 


s 


-3.085 


-4.090 


(2 


3) 


s 


+0.017 


-1.754 


(0 


2) 


A 


-6.767 


-5.719 


(1 


1) 


A 


-6.654 


-6.846 


(1 


2) 


A 


-2.427 


-4.692 


(1 


3) 


A 


-0.194 


-3.138 


(2 


0) 


A 


-7.172 


-7.356 


(2 


1) 


A 


-5.752 


-5.821 


(2 


2) 


A 


-4.457 


-4.169 


(2 


3) 


A 


-1.813 


-1.844 



